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First evidence for CP violation in two-body singly Cabibbo-suppressed decays of D° mesons reported by LHCb 
has recently aroused great interest in charm physics. In this document the latest measurements of CP violation 
in the charm sector are discussed. LHCb and CDF results on time-integrated CP asymmetries in D° — > n~ir + 
and D° -> K~K+ decays are presented in some detail. A search for CP violation performed by Belle in other 



two-body decays, namely D° — > K^n , D° 

Finally, results obtained by CDF with D° 
multi-body D decays, are shown. 



K° s 4'\ D+ 



decays, 



A7r+ and D+ 



r+i/M, is also presented. 



as well as by LHCb and BaBar with other 



1. Introduction 

The violation of CP symmetry is well established 
in the K° and B° meson systems [H-QI, an d first evi- 
dence in the system has been recently reported @ . 
Experimental measurements of CP violation (CPV) in 
the quark flavour sector performed so far are generally 
well described by the Cabibbo-Kobayashi-Maskawa 
mechanism d of the Standard Model (SM). How- 
ever, it is believed that the size of CPV in the SM is 
not sufficient to account for the asymmetry between 
matter and antimatter in the Universe hence ad- 
ditional sources of CP symmetry breaking are being 
searched for as manifestations of physics beyond the 
SM. 

The charm sector is a promising place to probe for 
the effects of physics beyond the SM. There has been 
a renaissance of interest in the past few years since ev- 
idence for D° mixing was first seen [ij [Io[ ■ Mixing is 
now well establishe d [ lll| at a level which is consistent 
with expectations jl2j . The recent evidence for CPV 
in singly Cabibbo-suppressed decays of D mesons to 
two-body final states, reported by LHCb [13|, has def- 
initely heightened the theoretical interest in charm 
physics. Prior to this measurement, CP asymmetries 
in these decays were expected to be very small in the 
SM [IllLzj], with naive predictions of up to 0(1CT 3 ). 
For this reason, the asymmetry measured by LHCb, 
characterized by a central value of 0(1O~ 2 ), came as 
big surprise. 

Unfortunately, precise theoretical predictions of 
CPV in this sector are very difficult to achieve, as 
the charm quark is too heavy for chiral perturbation 
and too light for heavy-quark effective theory to be 
applied reliably. This fact means we cannot conclude 
that the observed effect is a clear sign of physics be- 
yond the SM [l8rl2lj |. in order to investigate these 
promising hints further and clarify the picture, carry- 
ing out complementary measurements in other charm 
decays is of paramount importance. 

In the next sections, after a description of the LHCb 
and CDF measurements with two-body D° decays, 
we discuss other recent measurements by the BaBar, 
Belle, CDF and LHCb collaborations, involving two- 



body as well as three- and four-body D decays. 



2. D° 7r-7r+ and D° K K+ 

The time-dependent CP asymmetry Acp(f] t) for 
D° decays to a self-conjugate CP eigenstate / (with 
f — f) is defined as 



A CP (f; t) 



T(D°(t)^f)-T(D°(t)^f) 

r(D°(t)->/) + r(D°(t)->/)' 



where T is the decay rate for the process indicated. 
In general Acp(f; t) depends on /. For / = K~K + 
and / = 7r~7r + , Acp(f; t) can be expressed in terms 
of two contributions: a direct component associated 
with CPV in the decay amplitudes, and an indirect 
component associated with CPV in the mixing or in 
the interference between mixing and decay. 

The asymmetry Acp(f; t) may be written to first 
order as [22| 

A C p(f; t) = a d j P (f) + -a^, 

T 

where a^ip(f) is the direct CP asymmetry, r is the D° 
lifetime, and a^p is the indirect CP asymmetry, which 
is universal to a good approximation in the SM p3j . 
The time-integrated asymmetry measured by an ex- 
periment, Acp(f), depends upon the time-acceptance 
of that experiment. It can be written as 

A C p{f)=af P {f) + ^a 1 ^, 

where (t) is the average decay time in the recon- 
structed sample. Denoting by A the differences be- 
tween quantities for D° —5- K~K + and D° —5- n~Tc + 
it is then possible to write 



AA C P = A CP {K-K+) - A cp {it-it + ) 
= [a d c i P (K-K + ) a&(7r-7r+)] 



A(t) 



'■cp- 



In the limit of vanishing A(t), AAcp becomes equal 
to the difference in the direct CP asymmetry between 
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the two decays. However, if the time-acceptance is ^ 
different for the K~ K + and 7r~7r + final states, a con- 
tribution from indirect CPV remains. 

The LHCb collaboration performed a measurement 
of the difference in time-integrated CP asymmetries 



R-R+ and D° 



tt + using 



between D° 

0.62 fb^ 1 of data [13j]. The flavour of the initial state 
(D° or D°) is tagged by looking for the charge of the 
slow pion (7r+) in the decay chain D* + — > D°ir+. 

The raw asymmetry for tagged D° decays to a final 
state / is given by A raw (/), defined as 



-<4raw(/) 



N(D* + ^D (f)-Kf)-N(D*--,D {f)7v-) 
N(D*+^D°(f)Tr+)+N(D*-^-D°(f)n-)- 



where N(X) refers to the number of reconstructed 
events of decay X after background subtraction. 

To first order the raw asymmetries may be writ- 
ten as a sum of four components, due to physics and 
detector effects: 

Aaw(/) = Acp(f) + A D (f) + A d (tt+) + A P (D*+). 

Here, Ap>(f) is the asymmetry in efficiency for the D° 
decay into the final state /, Ad(7t+) is the asymmetry 
in efficiency for the slow pion from the D* + decay 
chain, and Ap(D* + ) is the production asymmetry for 
D* + mesons. The first-order expansion is valid since 
the individual asymmetries are small. 

For a two-body decay of a spin-0 particle to a self- 
conjugate final state there can be no D° detection 
asymmetry, i.e. Ap>(K~ K + ) = Apj(ir~Tr + ) = 0. 
Moreover, Ap,(tt+) and Ap(D* + ) are independent 
of / and thus those terms cancel in the difference 
A lavj (K~K + ) - A law (TT~ir + ), resulting in 

AA CP = A iaw (K-K+) - A iaw (TT-TT+). 

Note that the production asymmetry Ap(D* + ) can 
be neglected in the case of the CDF experiment. This 
is because, in contrast with LHCb which is a for- 
ward spectrometer and employs flavour-asymmetric 
pp collisions, CDF was a symmetric detector in pseu- 
dorapidity (77) and operated at a flavour-symmetric 
pp collider. Hence, integrating the measurement over 
a symmetric r\ range, the possible presence of a pro- 
duction asymmetry is removed by construction, i.e. 
A P {D*+) = 0. 

The mass difference (8m) spectra of candidates 
selected by LHCb as in Ref. [f3j], where 8m = 
m(hr /i + 7r+) — m(h~h + ) — m(n + ) for h = K, ir, are 
shown in Figure[TJ The D* + signal yields are approxi- 
mately f .44 x f 6 in the K~K + sample, and 0.38 x f 6 
in the 7r~7r + sample. 

Fits are performed on the samples in order to deter- 
mine A raw (K~ K + ) and A raw (7r _ 7r + ). The production 
and detection asymmetries can vary with px and pseu- 
dorapidity n, and so can the detection efficiency of the 
two different D° decays, in particular through the ef- 
fects of the particle identification requirements. For 
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Figure 1: LHCb fits to the 5m spectra, where the D° is 
reconstructed in the final states (a) K~ K + and (b) 7r~7r + . 



this reason, since the different masses of the K~K + 
and 7r~7r + final states lead to differences in the kine- 
matic distributions of accepted signal events in the 
two cases, LHCb performs the analysis in several kine- 
matic bins defined by the pt and n of the D* + candi- 
dates, the momentum of the slow pion, and the sign 
of p x of the slow pion at the D* + vertex. 

In each bin, one-dimensional unbinned maximum 
likelihood fits to the 8m spectra are performed. 
A value of AAcp is determined in each measure- 
ment bin as the difference between A raw (K~ K + ) and 
^rawfT - ^)- A weighted average is performed to 
yield the result AA C p = (-0.82 ± 0.21)%, where the 
uncertainty is statistical only. Including the system- 
atic uncertainty, LHCb measures the time-integrated 
difference in CP asymmetry between D° — > K~K + 
and D° —> tt~it~ 



Aj4pp Cb 



decays to be 
-0.82 ± 0.21(stat.) ± 0.11(syst.)] %. 



Dividing the central value by the sum in quadrature 
of the statistical and systematic uncertainties, the sig- 
nificance of the measured deviation from zero is 3.5o\ 
This is the first evidence for CP violation in the charm 
sector. 
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CDF Run II Preliminary 
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= [0.22 ± 0.24(stat.) ± 0.11(syst.)] %■ 
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Figure 2: CDF distributions of D°tt + mass with fit results 
overlaid for (a) D° -> 77-71+, (b) D° -> n~ir + , (c) D° -> 
K~K + and (d) L>° -> K~ K + . 



A similar measurement has been performed by 
CDF [2H ■ Figure [5] shows the invariant D°tt s mass 
for D° decays to 7r~7r + and K~K + , corresponding to 
an integrated luminosity of 9.7 fb -1 . The D* + signal 
yields are approximately 1.2 x 10 6 in the K~K + sam- 
ple, and 0.55 x 10 6 in the tt~tt + sample. These yields 
are very close to those of LHCb, despite a difference 
in integrated luminosity by a factor 15. 

CDF uses a different approach from LHCb to take 
into account differences in kinematics between the 
two decay modes. The kinematic distributions are 
equalized by means of an event- by- event reweighting 
technique, and then a single fit is performed on the 
reweighted sample integrated in phase space. The fi- 
nal result obtained by CDF is 



= [-0.62 ± 0.21(stat.) ± 0.10(syst.)] %, 



which deviates from zero by 2.7a. This result is com- 
patible with the LHCb measurement, with compara- 
ble accuracy and less than la difference between the 
central values. 

CDF has also provided individual measurements of 
A CP (D° ->• K~K+) and A CP (D° ^77-71+) [25], us- 
ing a carefully constructed combination of raw asym- 
metries measured from tagged D* + — > D°(K ~7r + )7r + , 
D*+ -> L>°(7r-7r+)7r+, D*+ -> D (K-K + )tt + and 
untagged D° — > K~ti + decays. The calculation as- 
sumes that the production asymmetry of D* + is negli- 
gible. This holds at CDF but not at LHCb, where this 
kind of measurement is considerably more involved 
and has not been performed yet. Using part of the 
full data sample, corresponding to about 6 fb -1 , CDF 
obtains: 



3. Other two-body decays 

The Belle collaboration has searched for CPV in the 
decay D° K%P° [26j[. where P° denotes a neutral 
pseudoscalar meson that is either a 7r°, tj, or if, using 
an integrated luminosity of 791 fb -1 . The observed 
KgP° final states are mixtures of D° — > K°P° and 
D° — > K°P° decays where the former are Cabibbo- 
favored and the latter are doubly Cabibbo-suppressed. 
SM K° — K° mixing leads to a small CP asymmetry 
in final states containing a neutral kaon, even if no 
CP violating phase exists in the charm decay. Fig- 
ure [3] shows the distributions of the mass difference 
M(D*) — M{D) for the various decay modes. No ev- 
idence for CPV in these decays is observed, as Belle 
measures: 



.4 
.4 
A 







^0 



■K"tt 



CP 
CP 

d°^k° sV ' 
cp - 



-0.28 ± 0.19(stat.) ± 0.10(syst.)] %, 
-0.54 ± 0.51(stat.) ± 0.16(syst.)] %, 
-0.98 ± 0.67(stat.) ± 0.14(syst.)] %. 



Belle has also searched for CPV in Df, -> (j)TT + [27j 
decays. For the </>7r + , this is achieved by measuring the 
CP violating asymmetries for the Cabibbo-suppressed 
decays D + — > K + K~tt + and the Cabibbo-favored de- 
cays D+ — > K + K~7T + in the K + K~ mass region of 
the (j> resonance, using 955 fb _1 of data. The mass 
distributions are shown in Figure |3J Belle finds about 
0.237 x 10 6 D ± and 0.723 x 10 6 Df decays. Assum- 
ing negligible CPV in Cabibbo-favoured decays, Belle 
measures 



.4 



CP 



= [+0.51 ± 0.28(stat.) ± 0.05(syst.)] %■ 



The result shows no evidence for CPV and agrees with 
SM predictions. 

Another relevant result has been obtained by Belle 
using 791 fb -1 of data, with the most sensitive search 
for CPV in the decays D + —> 7r + r/ and D + — > 
7T + r]' [28j |. Figure O shows the tt + t] and n+n' invariant 
mass distributions. The final results obtained by Belle 
for CPV in these decay modes are: 

Acp^ n+n = [+1.74±1.13(stat.)±0.19(syst.)]%, 



.4 



CP 



[-0.12 ± 1.12(stat.) ± 0.17(syst.)] %. 



A 



CP 



= [-0.24 ± 0.22(stat.) ± 0.09(syst.)] %, 



Again, no evidence for CPV is found with these 
modes. 



4. D° K°7r+7r- 

With the same trigger used to collect D° — > h + h~ 
decays, and using an offline selection based on a neural 
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Figure 3: Distributions of the mass difference M(D*) — 
M(D) for the D° -> KgP° decay modes studied by Belle. 
Left plots show the mass difference between D* + and D° 
and right plots show that between D*~ and D°. The top 
plots are for the Kgir final state, the middle plots for 
K%n, and the bottom plots for K%rf . The points with 
error bars are the data and the histograms show the results 
of the parameterizations of the data. 



network, CDF is able to reconstruct about 0.35 x 10 6 
decays to resonances of D* tagged decays, using 6 fb _1 
of integrated luminosity [2{|. Figure H] shows the 
KgTr + ir~ mass spectrum and the D° — > X^7r + 7r _ 
Dalitz plot obtained by CDF. The analysis is made 
both in a model-independent way, by binning the 
Dalitz plot and looking for bin-by-bin asymmetries, 
and by fitting the population of each resonance in D° 
and D° decays using an isobar model. 

The search for CPV in bin-by-bin asymmetries 
shows no evidence of deviation from a Gaussian dis- 
tribution with zero mean and unit width. The full 
Dalitz fit includes a parameterization of the efficiency 
over the Dalitz plane. An event- by- event reweighting 
is applied in order to equalize D* + and D*~ kinemat- 
ics. Again, no evidence of CPV is found in any of 
the considered resonant modes. CDF measures the 
following asymmetry integrated over the Dalitz plane: 

A C p = [-0.05 ± 0.57(stat.) ± 0.54(syst.)] %, 
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Figure 4: Invariant mass distributions with the fitted func- 
tions superimposed for (a) (jm + and (b) 4>ir~ , observed by 
Belle. 
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Figure 5: Invariant mass distributions for (left) n + rj and 
(right) n + n' final states, observed by Belle. Points with 
error bars and histograms correspond to the data and the 
fit, respectively. 



5. Other multi-body D decays 

LHCb has demonstrated the ability to select large 
samples of three-body decays. For example, about 
0.37 x 10 6 signal D+ R-K+TT+ decays with high 
purity have been selected using 35 pb _1 of inte- 
grated luminosity with data taken during 2010 [30] . 
The reconstructed K~K + n + mass distribution and 
the D+ -> K-K + n + Dalitz plot are shown in Fig- 
ure [7J The charge asymmetries in the control modes 
D + — > K~tt + tt + and Df — > K~K + tt + are investi- 
gated to eliminate the possibility of observing detector 
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CDF Run II preliminary, L = 6.0 fb 
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Figure 6: Invariant A"g7r + 7r" mass distribution (top) and 
D° -> Kstt + -r~ Daltiz plot (bottom) obtained by CDF. 
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Figure 7: Invariant K K + tt + mass distribution (top) and 
D + -> K'K +7T+ Daltiz plot (bottom) obtained by LHCb. 



asymmetries as signals of CPV. No significant charge 
asymmetries are observed, indicating that such sys- 
tematic effects are negligible at this level of precision. 
The Dalitz plot of the signal mode is then binned ac- 
cording to four binning schemes, two of which account 
for the resonant structure of the decay and two of 
which do not. The distributions of bin-by-bin asym- 
metries are consistent with no CPV in each of the 
binning schemes. 

With four-body decays the phase space becomes 
five-dimensional, and the difficulty of the analysis 
grows considerably. However, in four-body decays 
the method of T-odd moments becomes available. 
This has been applied e.g. by the BaBar collabora- 
tion searching for CPV in D + — > K^K + n + ii~ de- 
cays |3lj . The procedure consists in defining a triple 
product of the momenta of three of the particles, i.e. 
Ct = Pk+ ' {Pn+ x Pit-)- It is then possible to define 
the quantity At for D + decays (and its analogue At 
for D~ decays) as: 

. _ T(C T > o) - r(c T < 0) 
T ~ r(c T > o) + r(c T < o) ' 

Then, if the quantity At — \{At — At) differs from 
zero, this is a sign of CPV. No evidence for CPV is 



found, as BaBar obtains: 

At = [-1.20 ± 1.00(stat.) ± 0.46(syst.)] %. 

6. Conclusions 

We have summarized some of the latest results in- 
volving two-, three- and four-body D decays at BaBar, 
Belle, CDF and LHCb. In the past months, first evi- 
dence of CPV in charm decays has been first obtained 
by LHCb with a 3.5a significance, then followed by 
CDF with 2.7a. This has triggered a big effort by the 
theory community, in order to understand whether 
any effect of physics beyond the SM was manifest- 
ing itself or whether the result could be explained in 
the framework of the SM. Precise theoretical predic- 
tions in this sector are very difficult to achieve unfor- 
tunately, so it is not yet possible to draw firm conclu- 
sions. Nevertheless, if the result is confirmed at more 
than 5a with increased statistics by LHCb, one ex- 
pects CPV to show up in other charm decays as well, 
and the pattern of asymmetries will help theorists to 
decode the overall picture and determine whether it is 
truly a SM effect. For the moment, all the efforts look- 
ing for CPV in other D decays besides D° — > ir + ir~ 
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and D — > K + K have been frustrated. But the 
search has just started. 
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